Effective medium theory for disordered two-dimensional graphene 
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We develop an Effective Medium Theory to study the electrical transport properties of disordered 
graphene. The theory includes non-linear screening and exchange-correlation effects allowing us to 
consider experimentally relevant strengths of the Coulomb interaction. Assuming random Coulomb 
impurities, we calculate the electrical conductivity as a function of gate voltage describing quantita- 
tively the full cross-over from the fluctuations dominated regime around the Dirac point to the large 
doping regime at high gate voltages. We find that the conductivity at the Dirac point is strongly 
affected by exchange correlation effects. 
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I. INTRODUCTION 

Since the experimental realization of grapheneii" the 
problem of understanding the electronic conductivity of 
graphene has generated a large experimental and theo- 
retical effort, given its fundamental importance and its 
technological relevance. The unusual transport prop- 
erties of graphene^ arise mostly from its "Dirac spec- 
trum" - a linear, zero-gap dispersion relation of chiral 
Fermion carriers (electrons and holes), with the charge 
neutral "Dirac point" defined as the singular point (with 
vanishing density of states) where the electron and hole 
bands touch. In this work we present a theory that is 
able to describe at a qualitative and semi-quantitative 
level most of the current bulk transport experimental 
results on exfoliated graphene. The theory is designed 
specifically to handle the disorder induced density in- 
homogeneities present in graphene and dominating the 
physics at the Dirac point where the average carrier den- 
sity is zero. The theory is general: it can be applied 
to systems other than graphene (as formulated, it ap- 
plies readily to graphene bilayers); it does not depend di- 
rectly on the source of the spatial inhomogeneities; and it 
can be used to calculate other graphene transport prop- 
erties such as thermo-transport^ and magnetotransport. 
The theory takes into account non-linear screening ef- 
fects and exchange-correlation terms via a Hartree-Fock 
Local Density Approximation (LDA); in this sense the 
theory is non-perturbative and therefore, within the lim- 
its of this approximation, is valid for large values of the 
graphene's fine structure constant r s = e 2 /hvF£, where 
vf is the Fermi velocity and e the dielectric constant. 
For graphene charge transport, the theory is powerful 
being able to: 1) Predict values of conductivity at the 
Dirac point, a m in, in very good agreement with exper- 
iments; 2) Describe quantitatively the dependence of a 
on the doping cx V g including the crossover or plateau 
region, i.e. the region where V g is finite but smaller than 
the range of values for which cr(V g ) exhibits a linear be- 
havior; 3) Show that many-body terms are essential to 
understand the transport properties of graphene close to 
the Dirac point; 4) Explain the observed dependence of 



a m i n on the sample quality; 5) Explain the dependence 
of o~ m i n on r s observed experimentally- and predict its 
behavior over a wide range of r s , a prediction that could 
be tested experimentally; 6) Show that at large doping 
the spatial density fluctuations do not modify the lin- 
ear dependence of aonV g ; and 7) Predict the transport 
properties at very low but finite long-range disorder; a 
regime not accessible in previous theories. 

The success of the transport theory presented in this 
work relies on the proper treatment of the effects of the 
spatial carrier inhomogeneity. In a metal or semiconduc- 
tor, defects play a double role 6 : they act as scattering 
centers and locally modify the conduction-band carrier 
density. A defect effectively shifts the local Fermi level 
away from its average value, which has dramatic conse- 
quences on the properties of graphene close to the Dirac 
point where the average density, (n), is zero and there- 
fore the spatial density fluctuations completely dominate 
the physics. The source of the disorder— can be of differ- 
ent nature: short-range scatterers, such as atomic defects 
in the graphene lattice structure, ripples 8 or long range 
scatterers such as charged impurities, but all will induce 
spatial inhomogeneities of the carrier density that close 
to the Dirac point result in the formation of electron- 
hole, e-h, puddles^ ' 10 ! 11 ' 12 . This situation has emerged 
as the one realized in current experiments on exfoliated 
graphene. The ideal case of zero average density and zero 
amplitude of the density fluctuations has been considered 
in previous works that predicted a m in to be either or oo, 
depending on the nature of the disorder, or, in the limit of 
vanishing disorder, finite and universa l 13 i 14 i 15 ' 16 i 17 ] 18 i 19 . 

The formalism that we present does not depend on 
the source of the inhomogeneities, however recent exper- 
imental results ^ 20 ' 21 provide convincing evidence that 
random charged impurities - located in the graphene 
environment - are the dominant source of disorder in 
graphene. Therefore for concreteness we assume the car- 
riers in graphene to be subject to a disorder potential due 
to a 2D random distribution of uncorrelated Coulomb im- 
purities with surface density rii mp , placed at a distance 
d from the graphene layer. The paper is organized as 
follows: in section |TT] we develop and justify the effec- 
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tive medium theory for graphene in presence of a 2D 
random distribution of charged impurities, in section HTT1 
we present the results for the conductivity of disordered 
graphene and finally in section HVl we summarize the main 
results and discuss the range of validity of the approach 
presented. 



II. THEORETICAL MODEL 

As said in the introduction, in the remainder of the 
paper we only consider the disorder induced by random 
charged impurities in the graphene environment, moti- 
vated by the strong experimental evidenc o 20 ' 21 that this 
type of disorder dominates in current experiments on ex- 
foliated graphene on Si0 2 substrate. In order to be able 
to correctly describe the properties of graphene close 
to the Dirac point an accurate characterization of the 
strong, disorder induced, density fluctuations is neces- 
sary. This requires the calculation of the ground state 
for interacting massless Dirac electron in presence of dis- 
order, a very challenging problem. In absence of disor- 
der one can use a Kohn-Sham-Dirac Density-Functional- 
Theory, DFT, that appropriately takes into account the 
fact that in graphene the low energy quasiparticles be- 
have as chiral massless Dirac fermions^ 2 -. In this approach 
many-body effects enter via a local-density approxima- 
tion exchange-correlation potential. In presence of dis- 
order the DFT-LDA approach is computationally very 
expensive and as a consequence can only be used to cal- 
culate the ground state for few disorder realizations 22 
and does not allow for a statistical characterization of 
the carrier density fluctuations that requires the calcula- 
tion of the disorder averaged quantities. A computation- 
ally cheaper approach is the microscopic Thomas-Fermi- 
Dirac, TFD, theory that recently some of us have devel- 
oped and used to study the electronic structure of dis- 
ordered graphene^ 2 -. The TFD approach is very similar 
to the DFT-LDA one, the difference being that in the 
TFD theory also the kinetic energy is approximated by 
a local-density functional, whereas in the DFT-LDA the 
kinetic term is a differential operator that operates on the 
Kohn-Sham orbitals. For single disorder realizations the 
TFD results can be directly compared to the DFT-LDA 
results and we found them to be in very good agreement 
for impurity density values, rii mpi relevant for current 
experiments on exfoliated graphene: ni mp > 10 9 cm~ 2 . 
The Thomas-Fermi theory in general is accurate when 
the ground-state density varies on length scales bigger 
than the the Fermi wavelength, Af, i.e.: 
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yen inequality |T]) can be rewritten as: 
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Inequality ([2]) is equivalent to say that in any region with 
homogeneous carrier density the number of carriers has to 
be bigger than 1. As shown below the non- linear screen- 
ing properties of graphene guarantee that in presence of 
charged impurities this condition is satisfied for graphene 
also close to the Dirac point. 

To summarize the main properties of the graphene car- 
rier distribution close to the Dirac point in the lower 
color plot of Fig. Q] we show a typical result obtained us- 
ing the TFD theory for a given disorder realization with 
nimp — 5 x 10 11 cm~ 2 , d = 1 nm and r s = 0.8, values 
that are typical for current exfoliated graphene samples 
and accurately estimated by comparing the Boltzmann 
theory and the experimental results at high doping 20 . 
This plot exemplifies the main qualitative properties of 
the carrier distribution close to the Dirac point: I) the 
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FIG. 1: (Color online). n(r) lower plot, and corresponding 
<r(r), middle plot, at the Dirac point for a single disorder 
realization: rii mp = 5 x 10 11 cm" 2 , d — 1 nm r s — 0.8. 
Top panel shows a schematic of the Effective Medium used to 
describe the strongly inhomogeneous state of graphene close 
to the Dirac point. 



distribution is characterized by e-h puddles in quantita- 
tive agreement with surface probe experiment o 11 ' ; II) 
the typical size of an e-h puddle, defined as a region with 
same-sign charges, is of the order of the sample size, L, 
as expected for a semimetal close to the neutrality point; 

III) n(r) is characterized by two types of inhomogeneities, 
namely, wide regions (i.e. big puddles spanning the sys- 
tem size) of low density (an example is shown in Fig 1 
with the blue contour that contains « 10 electrons); and 
few narrow regions of high density (an example is shown 
in Fig 1 with the white contour containing 2 electrons); 

IV) the narrow puddle regions (peaks, dips) have a cor- 
relation length £ that depends on n lmp , r s and d^ For 
the parameter values used in Fig. [TJ typical for current 
experiments, £ « lOnm. The combination of the high 
density in the peaks/dips and the fact that in the low 
density regions n(r) varies over scales much bigger than 
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£ guarantees that the inequality ([2]) is satisfied over the 
majority of the graphene sample and therefore justifies 
the TFD theory. The properties I)TV) are essential for 
the development and justification of the effective medium 
theory, EMT, for graphene presented in this work. 

In presence of disorder, the problem of graphene trans- 
port then becomes one of how to properly take into ac- 
count the strong density inhomogeneities shown in Fig.[T] 
At high doping the graphene carrier density can be as- 
sumed to be homogeneous. In this regime, a has been 
shown experimentally to scale linearly with (n) . The the- 
ory that has been able to describe most successfully, away 
from the Dirac point, the large number of experimen- 
tal results for graphene transport is scattering from ran- 
dom charge impurities ^ 24 ' 25 ! 26 , where, within the RPA- 
Boltzmann approximation, we have: 
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where F(r s , d) was given in Ref. lid An important con- 
ceptual step to derive a correct transport theory valid 
close to the Dirac point is to introduce a local spatially 
varying "puddle" conductivity cr(r). The local conduc- 
tivity is a well defined quantity as long as <r(r) varies on 
length scales that are bigger than the mean free path I 
i.e.: 



Vcr(r) 



> i, 
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Because, as shown in Fig. Q] close to the Dirac point 
the density profile is characterized by large puddles with 
n 7^ the RPA-Boltzmann theory is expected to be valid 
locally inside the e-h puddles. To estimate a in a puddle 
we therefore use Eq. (|3|) replacing (n) by its local value 
n(r) so that in Eq. Q we can replace Vcr(r)/cr(r) with 
Vn(r)/n(r). Considering that / = ha/(2e 2 kp), using 
Eq. for disordered graphene inequality ^ takes the 
form: 
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imp, , sj . For r s = 0.8, d = 1 nm, 
11 cm -2 (parameter values relevant for current 



For graphene on Si0 2 (r s = 0.8) l/F{r s ,d) « 10. Close 
to the Dirac point there are few narrow regions of high 
density and size £(d, ni mp , r s ) 
"imp ~ 10 

experiments on exfoliated graphene) is £ ~ 10 nm< I as 
also shown in the lower panel of Fig. [TJ However, the re- 
gions where the carrier density changes over length scales 
of the order of £ are very sparse and the density landscape 
is characterized by wide regions where the density is ap- 
proximately constant on much bigger length scales, of the 
order of the system size. Because the small puddles of 
size rs £ are isolated (i.e. do not form a path spanning 
the whole sample) and occupy a small area fraction (we 
estimated it to be less than 20%i2) we can neglect their 
contribution to the graphene conductivity. In addition, 



given their high density, steep gradients at the bound- 
aries, and the fact that £ <C I, the small puddles should 
provide very small resistance to the carriers' motion. In 
the large regions where n(r) is almost uniform n is ap- 
proximately equal to the root mean square of the den- 
sity distribution n rms ~ rii m p 10 ^ 12 , and as a consequence 
in these regions the inequality ^ is satisfied for exper- 
iments on bulk exfoliated graphene for which we have 
L ~ l^iin and n imp is in the range [10 11 — 10 12 ] cm~ 2 . 
We conclude that the non-linear screening properties of 
graphene that justify the use of the TFD theory, and the 
linear dependence of a on n (Eq. ||3J)), allow for the def- 
inition of a local conductivity for graphene at the Dirac 
point. The conductivity spatial profile corresponding to 
the density profile of the lower color plot of Fig. []] is 
shown in the upper color plot of the same figure. 

Close to the Dirac point graphene can then be thought 
of as a highly heterogeneous material for which the EMT 
should be appropriate to calculate its transport proper- 
ties. Even though the EMT was proposed long ago for 
heterogeneous materials 6 -^, the development and justi- 
fication for graphene is not obvious and is the central 
result of this work. To derive the EMT one considers 
a single homogeneous region (e.g. area encircled by the 
white perimeter in the color plot for a in Fig. [1]) with 
conductivity a embedded in a uniform medium with ef- 
fective conductivity, oemt and then requires the spatially 
integrated inhomogeneities of the electric field E, due to 
the inclusion of the region with a ^ oemt, ^° ^ e ec l ua l 
to zero. For a 2D system this requirement translates into 
the equation: 
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The EMT is valid if the conductance across the regions 
where a is homogeneous is higher than the conductance 
inside these regions. For graphene at the Dirac point 
this requires that the conductance across the disorder 
induced p-n junctions, PNJ, (shown as black lines in the 
lower color-plot of Fig. [TJ) is higher than the conductance 
inside the puddles. Close to a PNJ Eq. ([3]) is not valid. 
If D is the length scale over which n changes across a 
PNJ and kp = y/im is the Fermi wave vector at the 
sides of the PNJ, then using TFD results for n(r) we find 
kpD ~ 1; for single disorder realizations this result can 
also be inferred from the DFT-LDA calculations^ 2 -. In 
this limit, the conductance per unit length of a PNJ,(?pnj 
is given b y 28 ' 29 : 
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Let p be the perimeter of a typical puddle and T its form 
factor then the condition that the conductance across the 
PNJ must be bigger than the conductance of a puddle 
takes the form: 
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where n is the density at the sides of the PNJ: n ~ 
n rms ~ ni m ? j 10 i 12 . Given the macroscopic size of the 
majority of the e-h puddles the length of the PNJ, p, 
is of the order of the sample size L ~ 1/xm, and then 
the inequality © is satisfied. We can therefore conclude 
that the EMT developed here properly describes the elec- 
tronic transport properties of bulk disordered graphcnc 
samples and that close to the Dirac point, transport in 
bulk graphene is dominated by scattering processes in- 
side the puddles and not across their boundaries as was 
assumed in previous works^S. Disorder averaging Eq. © 
we find: 
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where P(cr) is the probability for the local value of a. 
Using Eq. ([3|) we have: 
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a(n) - q-emt 
a(n) + oemt 



P[n] = 0. 



(11) 



To obtain an accurate value of ctemt it is crucial to 
select the correct form of P[n}. Choosing P[n] to be a 
Gaussian distribution of width n rms , for o~ m in we find 
ct emt ~ (e 2 /h)[l.U60/F(r s ,d)}[n rms /n imp \. Using a 
Lorentzian approximation for P[n] with width n^, we 
find cr^ MT = (e 2 /h)[2/F(r s ,d)](n L )/(n imp ). These ana- 
lytical results allow for a direct comparison with other 
results for a m in in the literature^. In our approach 
P[n] is accurately calculated using the TFD theory^ 2 - that 
takes into account nonlinear screening, exchange and cor- 
relation contributions 2 ^!, and can handle the crossover 
plateau region where P[n] has a bimodal character that 
is not captured by the Gaussian and Lorentzian approxi- 
mations. Only close to the Dirac point and for small val- 
ues of r s do the Gaussian and Lorentzian approximations 
give results that are comparable to the ones obtained us- 
ing P[n] calculated with the TFD theory, P[n] TFD . The 
TFD approach allows us to go to higher values of r s and 
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thereby also providing a range of validity for 



analytical theories. 



III. RESULTS 

In this section we present the results for the conduc- 
tivity of disordered graphene obtained using the effective 
medium theory developed in the previous section and 
the probability distribution, P[n] TFD , obtained using the 
TFD theory. Shown in Fig. [2] is a(V g ). We see that our 
theory correctly predicts a finite value of a at the Dirac 
point that is of the same order as that measured exper- 
imentally. At high gate voltages the linear scaling of a 
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FIG. 2: (Color online), a as function of V g for 3 different val- 
ues of nimp ■ The dashed (solid) lines are the results without 
(with) exchange. The dots at V g — are the results presented 
in Ref. [l(J for the same values of nimv- 



as a function of V g is recovered. One important achieve- 
ment of our theory is the ability to correctly describe the 
crossover of a from its minimum, <T m in, at V g = 0, to 
its linear behavior at high gate voltages. Moreover our 
work shows the importance at low gate voltages of the 
exchange term. In Fig. [3] we present the result for o~ m i n 
as function of ni mp . We see that our approach explains 
the non-universality of o~ m in from sample to sample and 
that exchange reduces the dependence of a m i n on n; TOp 
for impurity densities in the range [10 10 — 10 12 ] cm -2 . For 
nimp ;$ 10 10 cm -2 a m i n grows quite rapidly as ni mp is 
decreased. The same qualitative dependence of a m in on 
nimp has been recently observed experimentally as shown 



in Fig. 5 (a) of Ref. 1 2 ll. In Fig. [4] the dependence of a m i n 
on d is shown for a given value of ni mp . We see that a m i n 
depends weakly on d and this dependence is even weaker 
if exchange-correlation effects are taken into account. 
In Fig. [5] we show the results for a m in as a function of 
r s . The solid (dashed) line shows the values of a m in ob- 
tained including (neglecting) exchange. We see that cr m i„ 
has a non-monotonous behavior due to the fact that r s 
strongly affects both the carrier density spatial distribu- 
tion (by controlling the strength of the disorder poten- 
tial, screening and exchange); and the scattering time r. 
From the RPA-Boltzmann theory in presence of charged 
impurities we have that r = ^/n/(TrF(r s ,d)vFni m p), so 
that r depends on r s only through F(r s ). In Fig. [6] we 
show 1/F(r s ) for d = 1, to show how r decreases as r s 
increases. On the other hand as r s increases the width of 
P[n], i.e. the number of carriers, increases. For small val- 
ues of r s (< 0.1), as r s increases, the decrease of r dom- 
inates and Umin decreases very rapidly. In this regime 
the results without exchange are qualitatively similar to 
the results with exchange. For r s > 0.1 the effect of r s 
on the amplitude of the spatial density fluctuations be- 
comes very important. In Fig [7] the root mean square of 
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FIG. 3: (Color online). u min as a function of n imp , r s — 
0.8 and d — 1 nm. Solid (dashed) line are the results with 
(without) exchange. For comparison the results obtained in 
Ref. LLC] are also shown, dot-dashed line. 
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FIG. 5: (Color online). <r min as a function of r s for d = 1 
and rii mp = 10 11 cm -2 . 
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FIG. 6: 1/F(r 3 ,d = 1 nm) as a function of r s 



FIG. 4: (Color online). Dependence of <j m in on d for n; mp = 
10 11 cm -2 and r s — 0.8. Solid (dashed) line are the results 
with (without) exchange. 



the density fluctuations is plotted as a function of r s . We 
see that in presence of exchange n rms grows much more 
slowly than when exchange is neglected. This is due to 
the fact that in graphene, contrary to parabolic band 
materials, close to the Dirac point, exchange penalizes 
density fluctuations. As a consequence for r s > 0.1 the 
behavior of a m in is qualitatively different when exchange 
is taken into account. In Ref. [5J, r s was varied from 0.56 
to 0.8 and <J m in was found to remain constant in agree- 
ment with the results shown in Fig. [5] in presence of ex- 
change. We can then draw the important conclusion that 
to understand, even at a qualitative level, the transport 
properties of current graphene samples exchange contri- 
butions must be taken into account. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have derived and justified the effective 
medium theory for disordered bulk graphene. In particu- 
lar we have considered the case when the disorder is due 
to random charged impurities in the graphene environ- 
ment. Charge impurities have been shown to be the most 
important source of disorder for current experiments on 
exfoliated graphene. There are however other sources of 
disorder like ripple o 32 i 33 ' 34 i 35 and atomic defects 3 ^ that 
might play an important role in determining the trans- 
port properties close to the Dirac point of graphene sys- 
tems other than bulk exfoliated graphene. For example 
there is evidence that atomic defects might be the domi- 
nant source of disorder in epitaxial graphene. We should 
emphasize that our motivation was to calculate the in- 
trinsic conductivity of graphene close to the Dirac point. 
For this reason we have neglected the effects of contacts 
that have been shown to be import an t 37 ' 38 i 39 ' 40 but that 
do not modify the value of the graphene intrinsic con- 
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FIG. 7: (Color online), firms •) cLS cl function of r s for d = 1 
and riimp = 10 11 cm~ 2 . Solid (dashed) line are the results 
with (without) exchange. 



ductivity and in experiments can be minimized by per- 
forming four-terminal measurements. Our theory does 
not apply to graphene systems with a band gap^i and 
is only valid in the diffusive limit. To model the dis- 
order we have assumed the impurities to be randomly 
distributed in a 2D plane placed at a distance d from 
the graphene. The only two input parameters that enter 
the theory are the charge impurity density ni mp , and d. 
These two parameters have been accurately extracted by 
comparing the RPA-Boltzmann theory and the experi- 
mental results at high dopin g 20 i 21 . We have used the 
EMT developed here to study in details the transport 
properties of graphene close to the Dirac point, where 
the physics is dominated by the strong density fluctu- 
ations induced by the disorder. The parameters rii mp 
and d that enter the EMT should not be regarded as fit- 
ting parameters: to test our results for a given sample, 
rii m p and d should be extracted by comparing the value 
of the mobility at high doping with the one provided by 
the RPA-Boltzmann theory, and then used to obtain the 
EMT results at the Dirac point. A key ingredient for the 
EMT is the accurate calculation of the density probabil- 
ity distribution, P(n), close to the Dirac point taking into 
account non- linear screening and exchange-correlation ef- 
fects. We have used P(n) provided by the TFD theory^ 2 -. 
The TFD theory is justified for most of the current ex- 
periment on exfoliated graphene and gives results that, 
for single disorder realizations, are in a gree ment with the 
DFT-LDA approach developed in Ref.lsil On the other- 
hand due to its computational cost the DFT-LDA does 
not allow the calculation of disorder averaged quantities 
and therefore a statistical characterization of the density 
profile. The TFD-EMT approach that we have developed 
is valid when the inequalities (fTJ) , ([5]), © are satisfied. 
In these inequalities n is the density inside the electron- 
hole puddles. In the small regions, of size £ ~ 10 nm, 
because of the high carrier density, inequality (p} is al- 
ways satisfied. As we have shown, most of the graphene 



is covered by large puddles of size L p of the order of the 
system size, £, and with density of the order of n rms . As 
shown in Ref . [Iolfl2l n rms ~ rjriimp , where r\ is a number 
of order 1. Using these results we find that and © 
are equivalent apart from a factor of order 1 and that the 
inequalities ([1]), ((9]), take the form: 
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respectively. Given that for the cases of interest 
l/F(r s ,d) « 10, see Fig- [HI the second inequality in (TT5|) 
is the most stringent: the condition that the conduc- 
tance across the PNJ must be bigger than the conduc- 
tance inside the puddles, or equivalently the condition 
that the mean free path must be smaller than the pud- 
dles, sets the minimal size that the puddles must have for 
the EMT to be valid. Considering that L p ~ L > 1 /im 
and riimp S [10 10 — 10 12 ] cm~ 2 , we see that for bulk 
samples of exfoliated graphene the inequalities (fT2")l are 
satisfied. We conclude that for sufficiently disordered 
and large graphene samples the TFD-EMT approach is 
accurate to estimate the graphene transport properties, 
this conclusion is consistent with the results of Ref. |42| 
in which the graphene conductivity in presence of Gaus- 
sian disorder obtained using a full quantum mechanical 
calculation was found to be in agreement with the semi- 
classical Boltzmann one even at zero doping provided the 
disorder is strong enough. 

Using the TFD-EMT approach we find that at the 
Dirac point the theory gives a finite minimum con- 
ductivity in very good agreement with that observed 
experimentally 5,20,21 . In addition the theory is able to 
describe the crossover of the conductivity from its min- 
imum value at the Dirac point to its linear behavior at 
high doping. We have also calculated the dependence 
of Omin on riimp and found it to be in good agreement 
with recent experimental results 2 ^. We obtain the de- 
pendence of a m in on r s , finding that for r s < 0.1, a rn in 
decreases very rapidly with r s but for r s > 0.2, exchange- 
correlation terms are qualitatively important and must 
be included to correctly describe the minimum conductiv- 
ity. This result is in agreement with the experimental re- 
sults presented in Ref.[5|in which the minimum of conduc- 
tivity is found to be the same for r s — 0.5 and r s = 0.8. 
By comparing our results for a m in as a function of rii m p 
and r s to the experimental ones we find that exchange- 
correlation effects must be included when evaluating the 
transport properties of graphene at the Dirac point. The 
reason is that exchange-correlation effects, close to the 
Dirac point, at low riimp and high r s , strongly suppress 
the amplitude of the density fluctuations and therefore 
strongly affect the transport properties. The success of 
the EMT, coupled with the microscopic disorder-induced 
graphene electronic structure calculation, in obtaining 
the graphene conductivity near the Dirac point indicates 
that this technique should be useful in calculating many 
other properties of graphene in the theoretically difficult 
inhomogeneity-dominated regime near the charge neu- 
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trality point. 
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